Charge radius and dipole response of ^^Li 
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We investigate the consistency of the measured charge radius and dipole response of ^^Li within a 
three-body model. We show how these observables are related to the mean square distance between 
the ®Li core and the center of mass of the two valence neutrons. In this representation we find 
by considering the effect of smaller corrections that the discrepancy between the results of the two 
measurements is of the order of 1.5cr. We also investigate the sensitivity to the three-body structure 
of ^^Li and find that the charge radius measurement favors a model with a 50% s-wave component 
in the ground state of the two-neutron halo, whereas the dipole response is consistent with a smaller 
s-wave component of about 25% value. 
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INTRODUCTION 

The properties of the two- neutron halo nucleus ^^Li 
have been discussed in numerous theoretical and experi- 
mental papers but knowledge about its structure is still 
uncertain. In the past year, the results of two important 
measurements have been published, namely, the RMS 
(root-mean-square)charge radius of ^^Li obtained from 
laser spectroscopy [l| and the dipole response, which was 
probed by Coulomb dissociation on a Pb target |5| . The 
purpose of this paper is to investigate whether the two 
measurements can be explained simultaneously within a 
three-body model and to see what are the implications 
of the two results for the structure of ^^Li. 

The nucleus ^^Li is an excellent example of a so-called 
borromean system, a bound three-body system in which 
none of the two-body subsystems form a bound state. 
Thus ^^Li can be viewed as a three-body system consist- 
ing of a ^Li core and two valence neutrons where neither 
^°Li nor the dineutron system has a bound state. The 
nucleus ^^Li has only one bound state with a two-neutron 
separation energy of about 300 keV. 

Theoretical studies of ^^Li have primarily been based 
on three-body models of the two valence neutrons inter- 
acting with the ^Li core, and it has been a major chal- 
lenge over the past 10-15 years to obtain information 
about the neutron-core interaction, i. e., about the scat- 
teriiig states in the unbound nucleus ^°Li. Early stud- 
ies [3| assumed a dominant p-wave structure of the two- 
neutron halo, based on a rather high-lying p-wave res- 
onance in ^°Li. Another model Q assumed a shallow 
neutron-core potential, which does not have any bound 
states, and this resulted in a strongly s-wave dominated 
ground state of the two-neutron halo. In order to ex- 
plore the structure of ^^Li, a wider range of models were 
developed [H] and compared to measurements. 

A better calibration of three-body models for ^^Li be- 
came possible with an accurate measurement of the two- 



neutron separation energy, &„ = 295±15 keV [1], and 
a production measurement [7| which probed the contin- 
uum of ^"Li. The latter measurement suggested a p-wave 
resonance at about 540 keV and some influence of s-wave 
scattering near threshold. Let us also mention that the 
quadrupole moments of ^Li and ^^Li are the same within 
the 15% experimental uncertainty Q; this can be taken 
as a justification for using three-body models. 

An analysis of the /3-decay of ^^Li showed that 
about 45-55 % of the two- neutron halo must be in Pi/2 
orbits, whereas the remaining part would most likely oc- 
cupy s-waves. Analyses of the momentum distributions 
produced in high-energy breakup reactions also suggested 
a large s-wave component, from 20-40 % [10] to 35-55 % 
ll( . In fact, there seems to be a consensus toward a 
large s-wave component in two-neutron ground state, a 
component that is much larger than what was expected. 
This feature may be related to the famous parity inver- 
sion in the neighboring nucleus "'^^Be, where the ground 
state is a 1/2+ state and not a l/2~ state as one naively 
would expect for a p-shcU nucleus. It is of interest to see 
how the recent charge radius 1] and dipole response 
measurements fit into this trend. 



THREE-BODY MODEL INTERPRETATION 

There is a very close relationship within a three-body 
model between the charge radius and the dipole response 
of a two-neutron halo since they are both probes of the 
distance between the ^Li core and the center-of-mass of 
the dineutron system. The mean square charge radius 
{rp{Z,A)) for point-nucleons, for example, can be ex- 
pressed in terms of the charge radius of the core nucleus 
as follows, 



{rl(Z,A)) = {rl{Z,A-2)) + [jy (r^J, (1) 
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TABLE I: The measured charge radius of ®He [l^, the change 
in the mean square charge radius of ^^Li and ®Li, and the 
mean square distance between the *He core and the two- 
neutron halo in ®He are compared to results of three-body 
models (3BM) and GFMC calculations. The spin-orbit cor- 
rection, Eq. ([5]), was ignored and the charge radius of *He 
was set to 1.673(1) fm. 



Nucleus 


{r-Ly^' (fm) 


S{rL) (fm^) 


«2„) (fm^) 


"He exp. [15] 
3BM [14] 
3BM [17] 
GFMC [18] 


2.054(14) 
2.036 
2.011 
2.08(4) 


1.42(5) 

1.35 

1.25 
1.49(15) 


13.8 ± 0.5 

13.2 

12.3 
14.5 ± 1.5 



FIG. 1: Measured charge radii of He [TB] and Li P, [T^ iso- 
topes are compared to the predictions of the three-body mod- 
els (3BM) of '^He and "Li discussed in Ref. [H. 



where the second term is the correction which is caused 
by the center of mass motion of the core nucleus in the 
presence of the two valence neutrons. The correction is 
proportional to the mean square distance, {r^^2n)^ be- 
tween the core and the center of mass of the dineutron 
system. 

The total strength of the dipole response of a two- 
neutron halo nucleus is approximately given by the clus- 
ter sum rule 13], 



4 {rl2n), 



(2) 



which is also expressed in terms of the mean square dis- 
tance, {r'^^2n)^ between the core and the dineutron sys- 
tem. The sum rule assumes that the total dipole strength 
can be calculated by closure, which includes dipole transi- 
tions to the Pauli blocked core states. The effect of Pauli 
blocking is a minor but not insignificant correction, as 
we discuss in the next section. However, if we ignore it, 
we see that the charge radius and dipole response mea- 
surements are closely related, since they are both probes 
of the core-dineutron distance. 

To make contact with the measured charge radius rch 
we must correct the point-proton charge radii of Eq. 
([T|) for the finite sizes of protons and neutrons. We 
should also consider the so-called Darwin-Foldy contri- 
bution and the effect of the spin-orbit charge density dis- 



cussed in Ref. We therefore have the expression, 



3r 



4(mc) 



■ + {r%o, (3) 



where (i?^) = 0.757(14) fm^ and (Rl) = -0.1161(22) fm^ 
are the mean square charge radii of protons and neutrons 
[l^ . respectively, and 3?i^/[4(TOc)^] is the Darwin-Foldy 
term. 



Inserting Eq. ^ into Eq. ^ we obtain the following 
expression for the difference between the mean square 
charge radii of the two-neutron halo nucleus and the core 
nucleus. 



0.232 



(4) 



It is seen that the proton charge radius Rp in Eq. ([3]) 
drops out of Eq. ([4]), and so does the constant Darwin- 
Foldy term. The only two corrections that survive, 
— 0.232/Z and (r^)|° , are due to the non-zero, mean 
square charge radius of a neutron and the spin-orbit 
charge density of the two valence neutrons. 

In the following we ignore the spin-orbit correction to 
the charge radius except when otherwise explicitly stated. 
The reason is that this correction is model dependent, so 
it is not obvious how one should convert the measured 
isotope shift into a mean square, core-dineutron distance. 
One can calculate the spin-orbit correction in different 
models from the explicit expressions that are given in 
Ref. [l2| . In the shell model for spherical nuclei, with two 
valence neutrons occupying an unfilled {l,j) sub-shell, 
the spin-orbit correction is 



(5) 



where /i„ = -1.913 in the neutron magnetic moment, m 
is the neutron mass, and 

(l-s)=j(j + l)-^(^ + l)-3/4. 

In three-body models one would have to calculate (1 • s) 
numerically as an average value, since the 0^ ground 
state of the two-neutron halo contains many single- 
particle components [14i] . 

We show in Fig. [1] the measured RMS charge radii 
for the helium [isj and lithium [l|, [l^ isotopes. The re- 
sults for ^He and ^^Li are compared with the three-body 
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model calculations of Ref. 14|. These calculations em- 



ployed a density-dependent contact interaction to sim- 
ulate low-energy nn scattering, and the "^He-neutron 
Hamiltonian was calibrated to reproduce the known low- 
energy neutron-a scattering phase shifts. The calcula- 
tions for ^^Li are discussed below. 

The measured charge radius of ^He and the mean- 
square, core-dineutron distance we obtain when we ig- 
nore the spin-orbit correction in Eq. ^ are compared 
in Table I to the results of different models, namely, two 
three-body models of Refs. [ii,[i3|, and a recent GFMC 
(Greens function Monte Carlo) calculation [l^. The lat- 
ter is an improvement over the results that were pub- 
lished in Ref. [191] . 

Let us estimate the spin-orbit correction to the ^He- 
^He mean-square charge radius difference in the extreme 
limit where the valence neutrons occupy a pure (^3/2)^ 
configuration. Then (1 • s) = 1 , and we obtain from Eq. 
^ the correction (r-^)|° = -0.085 fm^. Another estimate 
is to evaluate the average spin-orbit correction in the 
three-body model developed in Ref. [l4|, since we know 
in this case the occupation probabilities of the single- 
particle orbits. The model quoted in line 5, Table II of 
Ref. 14], has 83% of the two-neutron halo in (^3/2)^ 
orbits. Considering all orbits of the halo we obtain the 
average value (1 • s) = 0.82. This implies the spin-orbit 
correction = -0.07 fm^, which would bring the 

GFMC calculation into perfect agreement with the mea- 
surement, whereas the three-body model [14] would be 
off by 10%, which is a discrepancy of almost 3a. 

The discrepancy between the measured and calculated 
charge radius of ^"'^Li which can be seen in Fig. [l]may re- 
flect uncertainties in the neutron-core Hamiltonian that 
we have used. The neutron halo ground state contains 
in this case 23% s-waves and is therefore referred to as 



the s23 model below. The model Hamiltonian [ij] was 
calibrated to reproduce the measured two-neutron sepa- 
ration energy [6| and also the p-wave resonance structure 
observed in Ref. 0, but there is still some uncertainty 
in the s-wave strength which we explore below. 



DIPOLE RESPONSE 

The measured dipole strength distribution [2] is com- 
pared in Fig. [2R to the prediction of an old three-body 
model of ^^Li [20|. The calculated distributions include 
the effect of the experimental energy resolution [2] . Al- 
though this model has a rather small two-neutron sepa- 
ration energy of 200 keV and did not include the recoil 
effects in a three-body system, it was able to reproduce 
fragmentation data at 800 MeV/nucleon fairly well 2l|, 
and Fig. [2]A_ shows that it also produces a dipole response 
that is in surprisingly good agreement with the measure- 
ment [3] . The strong peak near 300-400 keV (solid curve) 
is produced by the strong attractive interaction between 
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FIG. 2: The measured dipole response of ^^Li Q] is compared 
(A) to the calculations of Ref. i20j , with (solid) and without 
(dashed) the effect of the final state nn interaction, and (B) 
to calculations that are based on the s23 model described in 
the text. 



the neutrons in the final state. If this final state interac- 
tion is set to zero we obtain the dashed curve. 

The measured dipole strength distribution 0] is com- 
pared in Fig. ^fi to new calculations that are based on 
the s23 model [14|. The calculations (with and without 
the effect of the final state nn interaction) include the 
experimental energy resolution as done in Ref. [0]. The 
ground state of the s23 model has a realistic two-neutron 
separation of 295 keV and includes recoil effects in the 
three-body system exactly, c. f. the last term of Eq. (3.1) 



in Ref. [14| . The recoil effects are treated approximately 
in the three-body final state of the dipole response by 
ignoring the off-diagonal component pi ■p2/{Acm) (the 
last term in Eq. (3.3) of Ref. [3l), whereas the diago- 
nal term /(2Acto) (2Acm) is included through the 
reduced mass. With this approximation, the continuum 
dipole response can be computed with the method of Ref. 
[2 01 ] . We have checked the accuracy of this approximation 



with the discretized dipole strength function of Ref. [22] 
and have confirmed that the approximation works well 
for the ^^Li nucleus. The calculated peak (solid curve in 
Fig. [2J3) is higher and shifted slightly toward higher en- 
ergies in comparison to the data but the overall strength 
is very reasonable. 

The total dipole strength that was measured up to 
a 3 MeV relative energy is S(£;i)oxp = 1-42 ± 0.18 
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TABLE II: The difference between the measured mean square 
charge radii of ^^Li and ®Li [ij, and the mean square distance 
between the ^Li core and the two-neutron halo (extracted 
from Eq. for (r^)|'^=0), are compared to the results 

of three-body models (3BM), and the values extracted from 
the Coulomb dissociation (CD) experiment 0]. The assumed 
charge radius of ®Li was 2.216(35) fm. The last column shows 
the mean square distance between the two valence neutrons. 



Nucleus 




«2„) (fm^) 


{rl„) (fm^) 


"Li exp [1] 


1.175(124) 


37.9 ± 3.7 




revised [23] 


1.104(85) 


35.7 ± 2.6 




Old 3BM [3| 


0.728 


24.35 


39.0 


s05 3BM yj 


0.541 


18.7 


42.8 


s23 3BM yj 


0.789 


26.2 


45.9 


s32 3BM new 


0.895 


29.4 


51.6 


s50 3BM flO] 


1.120 


36.2 


70.1 


CD exp. [2] 


0.82(11) 


27.2 ± 3.5 





cluster rule. The result (Eq. (20) of Ref. Q) translates 
into an RMS core-dineutron distance of 4.73 fm. 

From an experimental point of view, extreme care must 
be taken when determining the total dipole strength asso- 
ciated with excitations of the halo, and when translating 
this strength into the size of the halo. This cannot be 
done accurately without some theoretical guidance be- 
cause the dipole strength can only be resolved at low 
excitation energies. Moreover, the cluster sum rule ([2]), 
which is sometimes used to determine the size of the halo, 
is not exact because it ignores the Pauli blocking of some 
of the final states as we discussed earlier. In the old 
three-body model of Refs. [3, [l^l, for example, the total 
strength obtained by numerically integrating the calcu- 
lated dipole strength distribution is 1.57 e^fm^, whereas 
the cluster sum rule strength is 1.73 e^fm^. That implies 
that the total strength is reduced by 10% compared to 
the cluster sum rule because of Pauli blocking. 



e^fm^. The calculated dipole strength up to 3 MeV, 
B{El,Erei < 3 MeV)cai, is 1.26 e^fm^ in the old model 
(Fig. EK) and 1.38 e^fm^ in the new s23 model (Fig. 
[2(3). Since the measured and calculated dipole strengths 
shown in Fig. [2] do not differ much, it seems reasonable 
to estimate the mean-square, core-dineutron distance as- 
sociated with the experiment by the simple scaling, 

~ B{El,Er.i<3MeYU ^"=■2"^'^^' 

in terms of the mean-square distance (?'c2n)3SM we ob- 
tain in the three-body model. We note that this scaling 
method is consistent with the cluster sum rule but it docs 
not necessarily require that the total dipole strength of 
the model is given by the cluster sum rule. We emphasize 
that the calculated dipole strength, which we insert into 
Eq. ([6]), is calculated in a model that respects the Pauli 
principle, whereas the cluster sum rule does not. 

The scaling method ^ gives essentially the same re- 
sult independent of which of the two models we use (the 
old model or the new s23 model). The average value is 
shown in the last line of Table II and it represents the 
mean-square, core-dineutron distance we extract from 
the measured dipole strength distribution. The last line 
of Table II also gives the difference in the mean-square 
charge radius of ^^Li and ^Li, which we derive from Eq. 
(HI) by ignoring the spin-orbit correction. 

We note that the RMS core-dineutron distance quoted 
in Ref. which is 5.01 ± 0.32 fm, is smaller than the 
5.22 fm value we obtain from the last line of Table II. 
The smaller size is the result of identifying the estimated 
total dipole strength (1.78 e^fm^) with the cluster sum 
rule, Eq. ([2]). An even smaller size was obtained in Ref. 
(2^ by fitting the measured dipole strength distribution 
of Ref. 0, and identifying the total strength with the 



CHARGE RADIUS MEASUREMENT 

_The measured charge radius of ^^Li is 2.467(37) fm 
1]. We emphasize that the uncertainty in the measured 
charge radius is partly due to the absolute calibration 
of one of the isotopes (^Li). The difference between the 
mean square charge radii for different isotopes is there- 
fore much more accurately determined. This is a great 
advantage for our discussion of the halo because the mean 
square distance between the ^Li core and the dineutron 
is directly related, according to Eq. to the differ- 

ence S{r'^f^) between the mean square charge radii of ^^Li 
and ^Li. We have considered this feature in our deter- 
mination of the uncertainties on the values of S{r'^i^) and 
(^c2n) shown in Table II. 

The change in the mean square charge radius from the 
reference nucleus ^Li was obtained from the measured 
isotope shift Sv'^g^ {A,7) and the calculated so-called fi- 
nite mass correction Si'^^{A, 7) according to the expres- 
sion [H 



"IS 



(7) 



1566.1 kHz 

The finite mass corrections that were used in Ref. [l[ have 



recently been reevaluated [23[ . Combining these new cor- 
rections with the measured isotope shifts of Ref. [l[ one 
obtains the 'revised' charge radius and core-dineutron 
distance shown in Table II. The important quantity to 
our discussion is the size of the halo which is here repre- 



sented by 



We see that the value we obtain from 



the Coulomb dissociation experiment is smaller than the 
values we obtain from the two interpretations [l], |2^ of 
the charge radius measurement. The deviation is in both 
cases a 2a discrepancy. 

In Table II we also give the results we obtain in var- 
ious three-body models of "'^^Li, ranging from the 'Old' 



5 



three-body model of Ref. Q to the s50 model of Ref. 
[lot , which gives a 50% s-wave component in the ground 
state of the two-neutron halo. It is seen that the s50 
model is in agreement with the charge radius measure- 
ment [li, whereas the s23 model of Ref. [l4| (with 23% 



s-waves in the halo ground state) is consistent with the 
CD experiment. 

Let us finally estimate the spin-orbit correction, Eq. 
dS]), which we have ignored so far when applying Eq. ([4]). 
In the extreme model, where the two valence neutrons 
occupy the configuration, the value of (1 ■ s) is 

-2, and from Eq. (P we obtain (r2)|° = 0.113 fm^. In 
the s23 model, where 61% of the halo is in the (pi/2)^ 
configuration, 23% are s-waves, and 16% are in higher 
{l,j) orbits, we obtain the average value (1 • s) = -1.09 
and {r^)2n = 0.062 fm^. Inserting this value into Eq. 
(HI, together with the core-dineutron distance obtained 
from the Coulomb dissociation experiment (last line of 
Table II), we now obtain the corrected value (5(r^;j)cD — 
0.88(11) fm^ for the difference between the mean square 
charge radius of ^^Li and ^Li. This implies that the 
charge radius of ^^Li extracted from the Coulomb dis- 
sociation experiment is 2.41(4) fm, which is consistent 
with the directly measured value of 2.467(37) fm [ij. 

The main part of the uncertainty in the charge radius 
of ^^Li stems from the uncertainty in the charge radius of 
^Li. A better representation of the discrepancy we obtain 
in our three-body-model interpretation of the measured 
charge radius and dipole response of ^^Li is the difference 



CD 



1.10(8) - 0.88(11) 0.22(14) fm^ 



(8) 



This is the result we obtain when we adopt revised finite 
mass corrections of Ref. [3] • The discrepancy is now of 
the order of 1.5cr. 



MATTER RADIUS 

Also quoted in Table II is the calculated mean square 
distance between the two halo neutrons, (j'^in), in ^^Li. 
This quantity is not probed by the two experiments dis- 
cussed above. It is seen that this distance increases dra- 
matically when the magnitude of the ground state s-wave 
component increases. 

The RMS matter radius, obtained from an analysis 
of interaction cross sections provides an additional con- 
straint on the size of the halo. The mean square matter 
radius of a two-neutron halo nucleus is determined by the 
size of the halo and the core nucleus as follows 
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FIG. 3: Error bands on the size of the two-neutron halo in 
^^Li obtained from the Coulomb dissociation experiment [3] 
(B(E1)), the charge radius 1] (rchargc), and the matter radius 
[24| (rmatter). The results obtained in three-body models, 
with 5 to 50% s-wave components, are indicated by diamonds. 



The RMS radii obtained in Ref. [24 are 2.43±0.02 and 
3.27±0.24 fm, respectively, for ^Li and ^^Li. From the 
halo distances given in Table II and the quoted matter 
radius of ®Li we obtain a ^^Li RMS matter radius of 3.29 
fm in the s23 model. Thus the s23 model agrees with the 
matter radius and also with the strength of the dipole re- 
sponse but the charge radius is too small compared with 
the measured value. The s50 model, on the other hand, 
agrees with the charge radius measurement, but the mat- 
ter radius (which is 3.66 fm) and the dipole strength are 
too large compared to experiments. 

The constraints on the mean square distances of the 
two-neutron halo obtained from the measurements of the 
charge radius, the dipole response, and the matter radius 
of ^^Li arc illustrated in Fig. [3] together with the pre- 
dictions of the three-body models. The limits from the 
charge radius are based on the revised values [23] in Table 
II and do not include the correction from the spin-orbit 
charge density. That correction may reduce the charge 
radius limits on (r^2n) by 1.9 fm^, according to Eq. ([l]). 



if we adopt the spin-orbit correction (r^)|° 
which we obtained in the s23 model. 

There is unfortunately some disagreement about the 
value of the matter radius which has been extracted from 
reaction data. An example is Ref. [2^ where RMS radii 
of 2.30 and 3.53±06 fm were obtained for ^Li and ^^Li, 
respectively. This constraint provides a lower limit which 
is close to the upper limit of the matter radius shown in 
Fig. 3, and it would therefore favor the s50 over the s23 
model. However, we do not think this is a reasonable so- 
lution because the s50 model produces a dipole response 
that has roughly a 33% larger strength than what has 



0.062 fm^ 
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been observed (compare the values of 
Table II). 
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We think that the l.Scr discrepancy we obtain in our 
three-body model interpretation of the measured charge 
radius and dipole response of ^^Li is most likely caused 
by the neglect of core polarization. Actually, it may seem 
surprising that the effect of core polarization is not much 
more dramatic. 

In order to estimate the effect of core polarization, we 
have performed Skyrme Hartree-Fock calculations for ^Li 
and ^^Li using the filling approximation and the SGII in- 
teraction. The results show that the mean-square charge 
radius of ^^Li increases by 0.3 fm^ from that of ^Li be- 
cause of the core polarization effect, which is caused by 
the proton-neutron interaction when the valence neu- 
trons occupy p-waves, whereas it increases by about 0.2 
fm'^ for s- waves. This accounts roughly for the dis- 
crepancy, Eq. ([8]), between the charge radius and the 
Coulomb dissociation experiment. 

We conclude that it would be desirable to extend the 
three-body model so that one can consider the effect of 
core polarization in a consistent way. Work in this direc- 
tion has already been done for the ground state of ^"'^Li 
[2^, and it is also being pursued by other groups 27 1. 



The work by Varga et al. [26j] shows that core polar- 
ization does play a significant role in their microscopic 
cluster model calculation of the charge radius of ^"'^Li. 
This can be seen in Fig. 2 of Ref. 1], where their re- 
sults, with and without the effect of core polarization, are 
compared to the measured charge radius. A further test 
of such models is provided by the measured quadrupole 
moments of ^Li and ^^Li and by the dipole response 
that was extracted from Coulomb dissociation data f^. 
In this connection, it would also be very useful to test 
the consistency of the measured charge radius and dipole 
response of ^Hc because the effect of core polarization 
should be much smaller for an a core. 

We are grateful to T. Nakamura for including the ex- 
perimental energy resolution in the calculated dipole re- 
sponse, and to S. Shlonio and S. C. Pieper for remind- 
ing us of the significance of the spin-orbit charge den- 
sity. This work was supported (H.E. and P.M.) by the 
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